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Critical velocity for vortex shedding in a Bose-Einstein condensate 
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We present measurements of the critical velocity for vortex shedding in a highly oblate Bose- 
Einstein condensate with a moving repulsive Gaussian laser beam. As a function of the barrier 
height Vo, the critical velocity Vc shows a dip structure having a minimum at Vo ~ /r, where is 
the chemical potential of the condensate. At fixed Vo ~ 7^, we observe that the ratio of Uc to the 
speed of sound Cs monotonically increases for decreasing cr/G where a is the beam width and ^ 
is the condensate healing length. We explain our results with the density reduction effect of the 
soft boundary of the Gaussian obstacle, based on the local Landau criterion for superfluidity. The 
measured value of Vc/Cs with our stiffest obstacle is about 0.4, which is in good agreement with 
theoretical predictions for a two-dimensional superflow past a circular cylinder. 
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I. INTRODUCTION 

A superfluid flows without friction but becomes dissi¬ 
pative above a certain critical velocity Vc via generating 
its elementary excitations such as phonons and vortices. 
The Landau criterion provides a conventional energetic 
consideration to determine the critical velocity, stating 
Vc = min[e(p)/p] [T], where e{p) is the energy of an ele¬ 
mentary excitation of momentum p. For a homogeneous 
system, the Landau critical velocity is equal to the speed 
of sound Cg. However, the dynamic response of a su¬ 
perfluid flow is significantly sensitive to the boundary 
condition of the system and hence a quantitative under¬ 
standing of the critical velocity has been a challenging 
task in the study of superfluidity. 

One of the paradigmatic situations considered in fluid 
mechanics is a two-dimensional (2D) flow past a circu¬ 
lar cylinder. For an incompressible flow, the local ve¬ 
locity is increased by a factor of 2 at the lateral sides 
of the cylinder [5] and the local Landau supersonic cri¬ 
terion suggests a critical velocity of Vc = 0.5cs that is 
independent of the radius R of the cylinder. Theoreti¬ 
cal studies showed that the onset of dissipation involves 
generating a counter-rotating vortex pair n a. More 
rigorous calculations, taking into account the compress¬ 
ibility of the superfluid and quantum pressure near the 
boundary of the cylinder, predicted that the critical ve¬ 
locity converges to Vc = 0.37cs in the large cylinder limit 
i? » C EHS], where ^ is the superfluid healing length. 
Experimental verification of the predictions on Vc/cg is 
highly desirable. 

In previous ultracold atom experiments, a similar situ¬ 
ation was investigated by stirring superfluid samples with 
a repulsive laser beam [T(1VHK] . The existence of finite 
critical velocities I14j and generation of vortex 

dipoles [n nails] were successfully demonstrated. How¬ 
ever, the measured values of Vc/cg ranged widely from 
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0.1 to 0.45, which did not allow a quantitative study of 
the homogeneous 2D problem. Theoretical investigations 
showed that the inhomogeneous density distribution of a 
trapped sample [TMIH] . 3D vortex dynamics [TM2T] . or 
the manner of stirring na na mHia should be criti¬ 
cal in the measurements. Recently, the vortex shedding 
dynamics was also investigated experimentally with po- 
lariton superfluids flowing past static defects [Ml ES] . 

In this paper, we systematically study the critical ve¬ 
locity for vortex shedding in a Bose-Einstein condensate 
with a repulsive Gaussian potential. We measure the 
critical velocity as a function of the barrier height Vq 
of the potential over a wide range of the beam width, 
10 < < 55. In particular, in order to address the 

2D homogeneous regime, we employ spatially large and 
highly oblate condensates, ensuring vortex dynamics in 
two dimensions. 

The key difference of a Gaussian potential from a hard 
cylinder is its soft boundary. A Gaussian potential, 
V{r) = Vo exp(—2r^/cr^), produces a density-depleted 
hole in the condensate when Vq > p, where p is the 
chemical potential of the condensate. The radius of the 
hole and the potential slope at the hole boundary are 
given as 


^HVo/p)/2 

(1) 

dV , 2 

— = ApR (7^, 

dr r^R 

(2) 


respectively. In comparison to the case with a hard cylin¬ 
der, the soft boundary reduces the density in the prox¬ 
imity of the obstacle and consequently lowers the local 
speed of sound. Then, from the local Landau criterion 
it is naturally suggested that the critical velocity of the 
hard cylinder defines an upper bound for that of the ob¬ 
stacle formed by the Gaussian potential. When the hole 
radius R becomes larger with higher Vq [Fig. 1(b)] and/or 
the beam width a decreases for fixed R [Fig. 1(c)], the 
obstacle would converge to the hard cylinder with stiff¬ 
ening its boundary. 

The main result of our measurements is that in the 
deep non-penetrable regime (i.e., Vq ^ p), the critical ve- 
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FIG. 1: (color online). Soft boundary of the optical ob¬ 
stacle formed by a repulsive Gaussian potential V{r) = 
Vb exp(—2r^/(T^). (a) Normalized potential slope at the ob¬ 
stacle boundary S = {j^/^)\dV/dr\r=R as a fnnction of a/^ 
and R/£,, where fj, and ^ is the chemical potential and the 
healing length of the condensate, respectively, and R is the 
obstacle radins such that V(R) — jj.. The boundary becomes 
stiffer as (b) R/^ becomes larger or (c) a/^ decreases. The 
corresponding trajectories are indicated by the arrows in (a). 


locity Vc increases with decreasing a/^ for fixed Vq/ fi and 
approaches about O.dcg. This observation is consistent 
with the expectation from the aforementioned discussion 
based on the local Landau criterion. Furthermore, the 
measured value of Vc/cs with our stiffest obstacle is in 
good agreement with theoretical predictions for a 2D su¬ 
perfluid flow past a hard cylinder. 
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FIG. 2: (color online). Vortex shedding in a highly oblate 
Bose-Einstein condensate, (a) Schematic of the experiment. 
An optical obstacle is formed by a repulsive Gaussian laser 
beam penetrating through the condensate and moves horizon¬ 
tally at a constant velocity in the center region, (b) An in-situ 
image of the condensate with the optical obstacle at the ini¬ 
tial position, (c) Probability P{v) of having vortex dipoles 
(blue circles) and the number of vortices (red diamonds) as a 
function of the velocity v of the optical obstacle. The barrier 
height Vo ~ 2.1fi and the Gaussian width a ~ 20^. The solid 
line is a sigmoidal function fit to P{v). The inset shows an 
image of a condensate containing a vortex dipole. Each data 
point was obtained from 15 realizations of the same experi¬ 
ment and the error bars indicate the standard deviations of 
the measurements. 


II. EXPERIMENT 

Our experiment starts with a Bose-Einstein condensate 
of ^^Na atoms in a harmonic trap formed by combin¬ 
ing optical and magnetic potentials m- The condensate 
fraction of the sample is over 90%. In a typical sam¬ 
ple condition, where the trapping frequencies are = 
27r X (9.0,400) Hz and the atom number of the conden¬ 
sate is Nq = 3.2(2) X 10®, the condensate healing length 
is ^ = h/y/2m^ ~ 0.46 /rm and the speed of sound is 
Cs — \/[ijm « 4.3 mm/s at the trap center, where H is 
the Planck constant divided by 27r and m is the atomic 
mass. By adjusting the trapping frequencies or the atom 
number of the condensate, ^ is varied up to 0.9 fim. The 
Thomas-Fermi radius and thickness of the condensate are 
Rtf/^ > 240 and ^tf/C < 6, respectively. In this highly 
oblate condensate, vortex line excitations are strongly 
suppressed [laisT] and the vortex dynamics is expected 
to be two dimensional. 

We adiabatically ramp up the power of a repulsive 
Gaussian laser beam in 1 s and hold it for 0.2 s to ensure 
that the condensate is stationary. Then, we translate 
the laser beam horizontally by 24 /im by using a piezo- 
driven mirror [Fig. 2(a)]. The velocity v of the laser beam 
is kept constant during the translation and controlled by 


adjusting the traveling time. The sweeping region of the 
laser beam is centered in the condensate [Fig 2(b)]. The 
density variation over the sweeping region is less than 
10% and the speed of sound can be well approximated 
to be spatially constant. After completing the sweeping, 
we slowly ramp down the laser beam power for 0.5 s, 
and take an absorption image of the condensate after 
expansion by releasing the trapping potential to detect 
vortices. 

The width a of the laser beam is calibrated from the in- 
situ images of very dilute samples penetrated through by 
the laser beam, taking into account the imaging resolu¬ 
tion of our system (« 5 /rm). The beam waist of the laser 
beam is 9.1(12) ^m at the focal plane and the Gaussian 
width a of the optical obstacle is controlled by defocusing 
the laser beam at the sample plane. The sample thick¬ 
ness (< 3 fj.m) is much shorter than the Rayleigh length 
of the laser beam and we ignore the beam divergence. 
The beam width is much smaller than the condensate 
radius and the chemical potential of the condensate is 
negligibly affected by the presence of the laser beam. 

A vortex dipole is identified with two density-depleted 
holes that are symmetrically located in the conden¬ 
sate with respect to the sweeping line of the laser 
beam [Fig. 2(c) inset]. After beinggenerated in the center 
region of the condensate, the vortex dipole moves toward 
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the edge of the condensate and splits into two individ¬ 
ual vortices that subsequently travel along the boundary 
of the condensate in the opposite direction. This pecu¬ 
liar obrit motion of the vortices was studied in Ref. m- 
A single vortex was observed occasionally at low v with 
probability < 4%. Because the vortex lifetime is over 
10 s [in], we attribute the single vortex to uncontrolled 
perturbations in sample preparation and we do not count 
it as a vortex dipole. 

The critical velocity Vc for vortex shedding is deter¬ 
mined from the probability distribution P{v) for having 
vortex dipoles after sweeping with the laser beam. Here 
P{v) is obtained from 15 realizations of the same exper¬ 
iment with a given sweeping velocity v, i.e., P(v) is the 
ratio of the number of images showing vortex dipoles to 
the total number of measurements. The critical veloc¬ 
ity Vc is determined by fitting a sigmoidal function to 
the probability distribution as P(v) = + 

[Fig. 2(c)]. We use the value of l.Sy as the measure¬ 
ment uncertainty of Vc, corresponding to the range of 
0.2 < P < 0.8. 


III. RESULTS AND DISCUSSION 

Figure 3(a) displays the results of the critical velocity 
as a function of the barrier height Vq for various beam 
widths. The critical velocity shows a dip structure hav¬ 
ing a minimum at Vq ^ fJ., clearly distinguishing the two 
regimes: a penetrable regime with Vq < fJ- and a non- 
penetrable regime with Vq > p.. The dip structure of the 
critical velocity can be accounted for by a consideration 
based on the local Landau criterion. For Vq < /r, the 
density minimum is located at the top of the Gaussian 
potential and decreases for higher Vq, and thus, lower¬ 
ing the local speed of sound. On the other hand, when 
Vb > /r, as described before, the potential slope becomes 
steeper with higher Vq and the density in the proximity of 
the obstacle boundary is gradually restored back to the 
bulk density, leading to a higher critical velocity. The 
dip structure around Vq = M becomes more pronounced 
with larger a/^ by lowering the minimum value [Fig. 3(a) 
inset]. 

In terms of the vortex shedding mechanism, the pen¬ 
etrable regime is different because there is no density- 
depleted region in the fluid and it has been anticipated 
that vortex nucleation would be initiated by generation 
of rarefaction pulses |8l|9l|28l|29]. Recently, the crit¬ 
ical velocity of penetrable obstacles for dissipation was 
investigated in effective ID systems |5DH55] . In one di¬ 
mension, the critical velocity is predicted to vanish as 
Vq approaches /x m but it is not the case in our 2D 
situation. Moreover, the critical velocity shows a quite 
intriguing dependence on a/^ [Fig. 3(a) inset]. Further 
investigation of the functional form of Vc{Vo /in 
the penetrable regime is warranted. However, in this 
work, we focus on the non-penetrable regime to address 
the 2D hard cylinder situation. 




FIG. 3: (color online). (a) Normalized critical velocity 
Vc/Cs versus the relative barrier height Vo/fi for various 
beam widths cr/G (c,0 = [24.5(24),0.46] (red diamonds), 
[16.0(14), 0.46] (black squares), [9.1(12), 0.46] (blue circles), 
and (9.1(12), 0.9] /xm (open squares). The inset displays Vc/cs 
at Vb « ^ as a function of cr/^ and the solid line is a line to 
guide the eye. (b) Same data in the non-penetrable regime 
(Vb > p) as a function of the hole radius R/^, together with 
additional data points for Vojfi > 10. 


In Fig. 3(b), we recast the data for the non-penetrable 
regime (Vq > /x) as a function of the hole radius R, to¬ 
gether with additional data obtained for Vq > lO/x with 
a = 20^ and 35^. Note that the hole rad ius is wea kly de¬ 
pendent on the barrier height as R (x ^). From 

the previous discussion of the boundary stiffness effect 
and S (X R/a"^ [Fig. 1(b)], one may expect that Vc/cs 
would be saturated to a certain value with increasing 
hole radius and also that the saturation behavior would 
be faster with smaller beam width. We see that the ex¬ 
periment data are roughly fit to the expectation in a small 
R region. However, when R is increased to be larger than 
CT, the growth rate of the critical velocity slows down and 
becomes negative. Note that at R/^ ^ 35, the critical ve¬ 
locity with small a= 20 is even lower than that with 
large cr/C = 54. 

One possible explanation for such climbing-over be¬ 
havior of the critical velocity is the imperfection of the 
laser beam profile. If the beam profile is not perfectly 
Gaussian, for example, the intensity profile of the outer 
part of the laser beam decays slower than exponential 
and then the potential slope S at the obstacle bound¬ 
ary would decrease with increasing R. We see that the 
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FIG. 4: (color online). Normalized critical velocity Wc/cs ver¬ 
sus a/^ for fixed VoZ/r ~ 7. Here a = 9.1(12) /rm for black 
closed squares and ^ = 0.46 /rm for black closed circles. The 
blue open circles show the results of theoretical calculation 
for a 2D hard cylinder of radius R = ay/]xi{yQjJi)j2 from 
Ref. [S]. The red open diamonds indicate the numerical re¬ 
sults for a Gaussian potential with Vq/m = 100 from Ref. [34] , 
where the potential slope S would be about two times higher 
than ours. The dashed line denotes the theoretically predicted 
value, Vc/Cs = 0.37 for a 2D homogeneous case with a hard 
cylinder in the large obstacle limit HE]. 

climbing-over of Vc occurs at i? ~ u in both of the data 
sets with cr/^ = 20 and 35. This seems to support the 
beam profile effect because in our experiment a is varied 
by defocusing the same laser beam. The factor of 
the laser beam is measured to be 1.2. 

To further investigate the soft boundary effect on the 
critical velocity, we take a different scanning trajectory in 
the parameter space of the Gaussian obstacle: decreas¬ 
ing a/^ for fixed Vq/m- In this setting, the hole radius 
R is also varied proportionally with a [Eq. (1)] and the 
trajectory corresponds to a diagonal line in Fig. 1(a). 
A scanning with fixed R/^, as depicted in Fig. 1(c), 
might be more ideal in terms of isolating the finite- 
R/^ effect OIHlIl], but this would require exponentially 
high Vo/fJ, for small a/^ as Vo/fj, = , necessar¬ 

ily recalling the outer part of the laser beam. We set 
Vq/^ « 7, where R ~ a and the normalized boundary 
slope S = {^/y)S « 4(cr/^)“^. In the previous measure¬ 
ments (Fig. 3), VciVo/fi) shows a maximum around this 
potential height. 

In the new set of measurements, we observe that the 
critical velocity monotonically increases as a/^ decreases 
(Fig. 4), which is consistent with our expectation from 
the boundary stiffness. It is worth noting that although 
the data are obtained from various samples with different 
healing lengths, they agree with each other in the plane 
of the dimensionless parameters a/^ and vdcs. This 
demonstrates the 2D character of the vortex dynamics 
in our system because the healing length is the only rel¬ 
evant length scale in 2D superfluid hydrodynamics. 


For our stiffest obstacle, Vc/cg is measured to be about 
0.40, which is slightly higher than the predicted value of 
0.37 for a 2D hard cylinder in the large-i?/^ limit. In our 
measurement, R > 10^ marginally satisfies the large ob¬ 
stacle condition, and the deviation of the measured value 
might be attributed to the finite-i?/^ effect. The depen¬ 
dence of the critical velocity on R/^ was investigated the¬ 
oretically EHin] and it was shown that Vc/cg gradually 
increases from the value of 0.37 as R/^ —>■ 0. In Fig. 4, for 
comparison, we display the theoretical results of Ref. [8] 
for a 2D hard cylinder of radius R = (T-\/ln(Vb/M)/2 (blue 
open circles) and the result of numerical simulations per¬ 
formed with a Gaussian potential with Vq/^l = 100 in 
Ref. |34j (red open diamonds). The experimental re¬ 
sults converge to the theoretical predictions when cr/^ 
decreases, i.e., the optical obstacle becomes similar to a 
hard cylinder by stiffening its boundary. 

Our results are inconsistent with the 1/i? dependence 
of Vc that was predicted from analytic analyses on the 
stability of a superfluid flow [BISS]- It is clearly seen 
that when R/d~ increases by a factor of 5, the 
Vc/cg decreases less than a factor of 2. Even without 
including the additional reduction effect due to the soft 
boundary of the obstacle, the critical velocity decreases 
much smaller than what would be expected from the 1/i? 
dependence. 

Finally, we want to recall a few aspects of the experi¬ 
mental condition that should be considered for the quan¬ 
titative comparison of the measurement results to theo¬ 
retical predictions for a 2D homogeneous case. First, the 
critical velocity was not measured for a steady flow con¬ 
dition but by sweeping a finite section of a trapped con¬ 
densate. The measurement can be affected by the sweep¬ 
ing manner. For example, if vortex nucleation requires 
a finite time, which might be longer than the sweeping 
time near the critical velocity 0112], it would result in 
a systematic, upward shift of the measured value of Vc- 
Second, although we expect suppression of 3D vortex dy¬ 
namics in a highly oblate condensate, it is an inevitable 
fact that the condensate has an inhomogeneous density 
distribution along the axial direction. Thus, one cannot 
completely ignore 3D responses of the condensate, in par¬ 
ticular, at the moment of vortex nucleation. It might be 
necessary or sufficient to introduce an effective speed of 
sound that would be lower than the peak value Cg of the 
condensate [5S] . 


IV. SUMMARY 

We have presented the measurements of the critical ve¬ 
locity for vortex shedding in highly oblate Bose-Einstein 
condensates and investigated the soft boundary effect of 
the moving obstacle formed by a Gaussian potential. Our 
results are consistent with a picture based on the local 
Landau criterion and the measure value of Vc/cg with the 
stiffest obstacle is in good agreement with the theoreti¬ 
cal predictions for a homogeneous 2D superflow past a 
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cylindrical object. This work has established a reliable 
experimental method to measure the critical velocity of 
a trapped condensate and its intriguing extension is to 
investigate the temperature dependence of the critical 
velocity, which might provide a new setting to study the 
role of thermal atoms in vortex nucleation [37H39]. 
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